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ABSTRACT 

We present the near- through mid-infrared flux contribution of thermally-pulsing asymptotic giant 
branch (TP-AGB) and massive red super giant (RSG) stars to the luminosities of the Large and 
Small Magellanic Clouds (LMC and SMC, respectively). Combined, the peak contribution from these 
cool evolved stars occurs at ~ 3 — 4 fj,m, where they produce 32% of the SMC light, and 25% of the 
LMC flux. The TP-AGB star contribution also peaks at 3 — 4 fim and amounts to 21% in both 
galaxies. The contribution from RSG stars peaks at shorter wavelengths, 2.2 /im, where they provide 
11% of the SMC flux, and 7% for the LMC. Both TP-AGB and RSG stars are short lived, and thus 
potentially impose a large stochastic scatter on the near-IR derived mass-to-light ratios of galaxies at 
rest- frame 1-4 /im. To minimize their impact on stellar mass estimates, one can use the M/L ratio 
at shorter wavelengths (e.g. at 0.8 - 1 /im). At longer wavelengths (>8 /im), emission from dust in 
the interstellar medium dominates the flux. In the LMC, which shows strong PAH emission at 8 /xm, 
TP-AGB and RSG contribute less than 4% of the 8 /xm flux. However, 19% of the SMC 8 /xm flux 
is from evolved stars, nearly half of which is produced by the rarest, dustiest, carbon-rich TP-AGB 
stars. Thus, star formation rates of galaxies, based on an 8 /im flux (e.g. observed-frame 24 fim at 
z = 2), may be biased modestly high, especially for galaxies with little PAH emission. 



Subject headings: galaxies: stellar content 
parameters 

1. INTRODUCTION 

Stellar masses of galaxies are typically estimated from 
a mo del mass-to-light (M/L) ratio (e.g. iBell fc de Jond 
1200 It) , and a measurement of the integrated luminos- 
ity at optical or near-infrared wavelengths. Until re- 
cently, galaxy M/L ratios were thought to be rela- 
tively well behaved, especially at near-infrared wave- 
lengths where the effects of dust and massive main se- 
quence stars are minimized. However, in the past decade 
this assumption has been shown to be flawed because 
of the influence of short-lived but extremely luminous 
thermally- pulsing asymptot i c giant branch (TP -AGB) 
stars (M araston et al.l 120061 : iZibetti et akl I2009L 0.8 < 
M < SMpi) and ma s sive red super giant s tars (RSG; 
iDalcanton et al.ll2012t iMelbourne et al.ll2012D . 

Just as the near-IR fluxes of galaxies are often used 
to estimate their stellar masses, mid-IR fluxes can be 
used to estimate their star formation rates. For instance, 
at rest-frame 8 /xm, the flux from starbursts is typically 
dominated by polycyclic aromatic hydrocarbon (PAH) 
emis sion which scales with th e star formation rate (SFR; 
e.g.. lDiaz-Santos et al.|[2008l ). Likewise the rest- frame 24 
/im flux has been shown to scale with the IR luminos- 
ity of a galaxy and is therefore used as a SFR indicator 
(|Ch arv fc Elbad I2OOII ). However, calibrations of these 
SFR indicators typically ignore the possibility of stellar 
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contamination in the flux at these wavelengths. Sev- 
eral works have shown that, under certain assumptions, 
AGB stars could contribute significantly to the mid-IR 
fluxes of galaxies and therefore affect the SFRs estimated 
from their mid-infrared fluxes (e.g.. iVerlev et al.l 120091 : 
Kelson fc HoidCTl [20Tot iChisari fc KelsonI I2012D . How- 
ever, as will be discussed in this paper, some of their 
assumptions may be flawed. 

The AGB is a short-lived phase (~ 1 Myr, e.g., 
IGirardi et all l20Tol) near the end of stellar evolution 
when a star undergoes helium shell burning and swells 
to a large size, producing a relatively cool photosphere 
and, potentially, a very infrared- (IR-) luminous star. 
Because of the cool temperatures, complex molecules 
and dust can form in the atmosphere f^e.g.. | Bowe: 



T988H\yinters et al ll200a[i553l:lSchirrmacher et al"ll200 
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Woitke|[200a [20071 : iMattsson et al|[2008l: Ivan Loon et all 
2008), which in turn can drive strong stellar winds (e.g., 
Sedlmayr & Dominik 1995; Elitzur & Ivczic 2001, and 
references therein ) that ultimately return as much as 
half of the stellar mass to the interstellar medium. An 
isolated star will start its AGB phase with an oxygen- 
rich atmosphere (e.g.,0-AGB). However, during the later 
stages of AGB evolution, a star undergoes a series of 
thermal-pulses and carbon, dredged up from the inte- 
rior, pollutes the atmosphere (Ibcn 1983). When the 
C/0 ratio exceeds unity the star is classified as a car- 
bon star (e.g., C-AGB), and complex molecules and dust 
grains form in the escaping winds. These dust grains 
tend to obscure the star at optical and near-IR wave- 
lengths, but the same dust emits in the mid-IR. Only 
20% of the AGB stars in the LMC/SMC are estimated to 
be thermally-pulsing (Bruzual, private communication). 
However, these TP-AGB stars can be very luminous in 
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the infrared, typically exceeding the luminosity of the tip 
of the red giant branch (TRGB), sometimes by several 
magnitudes. 

Red supergiants are helium-burning massive stars that 
form a tight sequence in optical and near-IR CMDs (e.g. 
Figured]). Stars on this seq uence are also termed red 
helium burning stars (RHeB | Dohm-Palmer fc Skillmaiil 
120021: iMcQuinn et all 120 lit iDalcanton et al.ll201lD . At 
the luminous end, RSGs are truly massive, 9 solar masses 
or more and thus were born in very recent star formation 
episodes (e.g. < 0.1 Gyr). At the luminosity of the K- 
band TRGB, the Padova isochr ones p redict more modest 
masses, 6-7 (iBertelh et al.ll200llMarigo et al.ll2008l: 
iBertelh et all [2009: Gir ardi et aLlFoiOD . and in fact the 
RHeB sequence continues to lower luminosities for stars 
of even smaller mass (e.g. IMcQuinn etHIIMTI ). 

Combined, TP-AGB and RSG stars have been shown 
to account for as much 40% of the integrated 1.6 fim flux, 
even in local galax ies, where the red giant branch (RGB) 
is well developed (|Melbourne et al.ll2012|) . In the early 
universe, when the RGB has not had time to form, the 
TP-AGB and RSG stars are expected to contribute even 
larger fractions (as much as 70%) of the near-IR light, 
even tho ugh they represent a t iny fraction of the ste l- 
lar mass (|Maraston et al.ll2006t iMelbourne et al.ll2012[ ). 
Thus, these stars must be accounted for to accurately 
model the infrared M/L ratios of galaxies (|Zibetti et al.l 
[2009; Ilbcrt et al. 2010). 

Not only do these stars contribute to the near- 
IR luminosities of galaxies, but they a lso can be 
significant sources of mid-IR light (e.g., iBover et al.l 
I2OII1) . TP-AGB stars, in particular, can produce 
significant amounts of w arm dust (e.g., iBowe 3 IT988I: 
iSchirrmacher et al]l2003[ ). These dusty shells both emit 
in the mid-infrared and obs cure the th e TP-AGB stars 
at optical wavelengths (e.g.. lReidlll991| ). 

The flux contribution of TP-AGB and RSG stars 
has now been well quantified at 0. 8 and 1.6 /xin in a 
set of 23 local dwarfs and spirals (|Girardi et al.l I20I0I : 
IMelbourne et al.ll2012l ). Similarly the 3.6 and 4.5 /xm 
fiux contribution of the most extreme (dusty) TP-AGB 
stars has been quantified in a set of six nearby galaxies 
(|Gerke fc Kochanekll20ll ). Likewise the mid-IR (8 - 24 
/im) fiux contribution from carbon stars has been esti- 
mated fro m stellar population syn thesis models of SINGS 
galaxies (|Chisari fc KelsonI 120121 ). However, there has 
not been a study that quantifies the fiux contributions of 
these stars across the full 1 — 24 ^m wavelength range for 
a galaxy where complete samples of TP-AGB and RSG 
stars are actually resolved and individually measured. 

In this paper, we use complete samples of TP-AGB and 
massive RSG stars in the Magellanic Clouds from the 
Spitzer legSiCy progra m "Surveying the A gents of Galaxy 
Evolu tion" (SAGE iBlum et al.i I2OO6I: iMeixner et al.l 



12051: IGordon et all [20111: IBover et alTl201lD . We build 
off of the results of IBover et al.l (j20lH) . who showed that 
the TP-AGB and RSG stars contribute significantly to 
the near- and mid-IR point-source fiuxes in these galax- 
ies. We now show that their contribution to the total 
(stellar and interstellar) fluxes of the SMC and LMC 
varies significantly from 1 to 24 /zm. We also find signifi- 
cant differences between our measured flux contributions 
and those from several recent studies in the literature. 
We find a significantly larger TP-AGB contribution at 



3.6 iim than is foun d in iGerke fc Kocharie3 (j2012f ) and 
iMeidt et al.l (|2012aD . Conversely, we fi nd the 8-24 nm 
AGB flux c ontributions estimated i n [Kelson fc Holdeiil 
(|2010f l and iChisari fc KelsonI (|2012l ) are likely signifi- 
cantly over-estimated. We discuss the implications of 
these results for estimates of stellar mass and SFRs of 
galaxies. 

2. THE DATA 

This paper uses several types of observations: (1) 
ground-based, seeing-limited near-infrared imaging to 
measure the near-infrared fluxes of TP-AGB and RSG 
stars in the MageUanic Clouds, (2) space-based, low- 
resolution near-infrared imaging to determine the total 
integrated near-infrared luminosities of the Magellanic 
Clouds, and (3) space-based, mid-infrared imaging to de- 
termine the mid-infrared fluxes of both the stars and the 
integrated galaxy light. In addition to the IR data, we 
use ground-based optical data to explore the flux con- 
tribution of short-lived massive main sequence stars at 
shorter wavelengths. Details of these datasets are pro- 
vided below. All magnitudes are reported in the Vega 
system. 

2.1. Stellar Photometry 

Multi- wavelength stellar photometry of the stars in the 
Large and Small Magellanic Clouds (LMC and SMC, 
respectively) was compiled by the SAGE team. J 
and Kg (ie., 1.2 and 2.2 /xm) near-infrared photome- 
try was o btained from the Two Micron All Sky Survey 
(2MASS: ISkrutskie et aLll2006D . Mid-infrared photom- 
etry at 3.6, 4.5, 5.8, and 8 /^m was obtained with the 
Spitzer Space Telescope InfraRed Array Camera (IRAC; 
Fazio et al. 200^). Photometry at 24 /im was obtained 
with the Mul ti-band Im aging Photometer for Spitzer 
fMIPS: iRieke et al.|[200l . A discussion of the Spitzer 
observations, image reductions, p hotometry, and sourc e 
matching for the SMC is given in IGordon et aP ([Ml 
and the LMC photometry is discussed in IMeixner et all 
f2006). Optical, BVI, data was compiled by the Magel- 
lanic Clouds Photometric Survey (MCPS; Zaritskv et al.l 
[2001 . 

Stellar clas sifications fo r stars in the SAGE data were 
presented in IBover et al.l (|2011f ). Here we briefiy sum- 
marize the method; for a full description of the stellar 
classification and sources of possible contamination see 
IBover et al.l (|2011D . Classification was carried out by ex- 
amining the locations of stars in multi-wavelength color- 
magnitude space. Less-dust-obscured TP-AGB and RSG 
stars were identified in near-infrared color magnitude di- 
agrams (CMDs), where they are among the most lumi- 
nous stars, e.g. above the i^s-band TRGB. At these 
wavelengths, sequences of carbon-rich and oxygen-rich 
AGB stars (C-AGB and 0-AGB stars respectively) sep- 
arate by J — Ks color, with C-AGB stars typically redder 
than 0-AGB stars. Likewise, RSG stars form a tight se- 
quence slightly blueward of the 0-AGB stars. These se- 
quences are clearly visible in the J — CMD of the 
LMC shown in Figure [T] which also shows the color- 
magnitude boundaries used to ident ify these different se- 
quences (adopted from iCioni et al.l [20061. While these 
boundaries are useful for identifying different classes of 
stars, we note that there will be some overlap between the 
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Fig. 1. — J — -R " s vs. Kcj CMDs for the LMC (grey density plots). Left panel overlays the Padova isochrones ijBertelli et al.l [20081 : 
IMarigo et al.|[2"00SI : IBertelli et ani2009IMGira rdi ct al. 2010). Center panel overlays a statistical es timate of the Milky Way foreground stars 
in the direction of the LMC from TRILEGAL (Girardi et al.ll200g) . The right panel shows the'Bover et al.' fSoTT) classification of RSG 
(yellow), O-AGB (blue), and C-AGB (red) stars, with selection regions defined by Cioni et al. (2006). A luminosity cutoff was made at the 
TRGB , below which the fluxes of these stars become negligible compared with the RGB . The dust enshrouded x- AGB (green) were selected 
by Boyer et al. (2011) at longer wavelengths; J — Ks > 3.1 mag or [3.6] — [8] > 0.8 mag and [8] < 10 mag. Because of their dusty envelopes, 
the x-AGB stars can be fainter than the ii'g-band TRGB and yet be significantly brighter than the TRGB at longer wavelengths. 



populations. Boyer et al. (2011) provides a detailed dis- 
cussion of classifications including estimates of the num- 
bers of misclassified stars, which are small. We note that 
the an omalous 0-rich AGB stars reported bv lBover et al.l 
are included with the O-AGB sample here. 

The most dust enshrouded TP- AGB stars are signif- 
icantly extinguished even in the near-IR, and the refore 
they m ust be identi fi ed at longer wavelengths (e.g. JReidl 
119911 ). iBover et al] H2OO90 found that short wavelength 
searches can miss up to 40% of the TP- AGB stars. These 
so-called extreme AGB stars (x-AGB) were identified by 
their colors at mid-IR wavel engths (3.(3 to 8 ^ m) . Most x- 
AGB stars are C-rich (van Loon et al.ir2005( ). Combining 
these two identification methods allows for the selection 
of near-complete samples of TP-AGB stars (e.g. these 
classifications do not miss the most dust enshrouded 
stars, which can be missed by classifications that rely 
entirely on shorter wavelengths). 

While these techniques will recover the bulk of the 
luminous TP-AGB stars, some TP-AGB stars will be 
missed because of they are at a minimum in their thermal 
pulsation cycle and are fainter than the TRGB. Models 
of the TP-AGB stars in the Magellanic Clouds based on 
the Padova stellar evolution tracks reveal that < 10% 
of them li e below the TRGB due to the thermal pulse 
cycle (e.g. IMarigo fc Girardill2007D . We have included 
these sub-luminous TP-AGB stars in the uncertainties 
of the integrated TP-AGB fluxes by assuming they have 
the luminosity distribution at the top magnitude of the 
RGB. We added this uncertainty in quadrature to the 
other sources of uncertainty, increasing the upper error 
bars in Table 1 and Figure 2 by < 3%. 

We also note that there are AGB stars that are not 
thermally-pulsing. These are lower-luminosity AGB 
stars and also lie below the TRGB. These stars cannot 
be easily selected in CMD space because they overlap 
the RGB stars. However, as these AGB stars are also 
relatively rare by number compared with the equally lu- 
minous RGB, they will not contribute significantly to the 
total flux of the galaxy. 

The total numbers of SMC and LMC stars within each 



stellar classification bin are presented in Table [TJ 

2.2. Foreground Contamination 

The CMDs of the LMC and SMC also contain some 
Milky Way foreground stars with the colors and lumi- 
nosities of Magellanic Cloud AGB and RSG stars. We 
applied a statistical correction for foreground contami- 
nation of Milk y Way stars at eac h wavelength. We used 
TRILEGAL ( Girardi et al.ll2005[ ) to model the numbers 
and fluxes of Milky Way stars in the direction of the Mag- 
ellanic Clouds in the same filters as our data. The TRI- 
LEGAL photometry was then run through the same clas- 
sification scheme as the RSG and TP-AGB stars. The 
contamination is largest for the RSG population of the 
SMC ('- 25%). In the LMC the RSG contamination 
is negligible because the higher metallicity of the LMC 
causes all of the RSG and AGB stars to be redder than 
they are in the SMC, resulting in less overlap with the 
foreground sources in the near-IR CMD (see Figure [1]). 
Likewise, the AGB populations in both galaxies are ef- 
fectively foreground contamination free. The estimated 
numbers of foreground stars within each stellar class are 
also presented in Table [H 

Note: stars in the Milky Way foreground clusters 47 
Tuc and NGC 362, in the direction of the S MC, were 
removed from the data (see IBover et al.l[20ll . 

2.3. Integrated Galaxy Photometry 

Measuring the total luminosity of the Magellanic 
Clouds is challenging because of their very large areas on 
the sky. Ground-based, near-IR images are dominated by 
thermal background fiux, and thus accurately discerning 
sky from unresolved galaxy becomes nearly impossible. 
Fortunately these issues are largely eliminated by going 
to space, where the IR thermal background is minimal 
and COBE/DIRBE (Hauser et al. 1998) provides a large 
beam on the sky (nearly a degree). 

We use the C0I3E/piRBE integrated fiux density es- 
timates of llsrael et al.l ()2010D for the Magellanic clouds 
at 1.25, 2.2, and 3.5 fim. The SMC fiux is measured 
within a 12?3xl2?3 area, and the LMC flux covers a 
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15?3x 15?3 area. These measurements were carefully cor- 
rec ted for fore g round (and background) contamination 
by [Israel et al.l (|2010[ ). The Israel et al. flux density es- 
timates for the Clouds are presented here in Table [1] 

Similarly, the global SMC and LMC IR fluxes from the 
Spitzer data at 3.6, 4.5, 5.8, 8, and 24 /im were made by 
[Gordon et al.l (|2011[ ) and Gordon et al., (in preparation). 
These fluxes were measured within circular apertures of 
radius 2? 25 for the SMC and 3? 75 for the LMC. For both 
galaxies, a sky annulus with a 0?2 width was used to 
estimate and subtract the background/foreground. To 
within the uncertainties, the integrated DIRBE fluxes 
at 3.5 /xm from Israel et al. match the Spitzer fluxes 
at 3.6 fim from Gordon et al. The Spitzer fluxes are 
provided in Table [1] 

The spectral energy distributions of the SMC and LMC 
are shown in Fi gure [2l Also ov erlaid are dust model 
predictions from lBot et al.l ()2010f ). While the stellar flux 
is declining with wavelength, the warm dust and PAH 
components are increasing. 

We compare these integrated galaxy fluxes gathered 
from the literature to the fluxes from the TP-AGB and 
RSG stars as measured from 2MASS and Spitzer pho- 
tometry. The 2MASS J and ifs-bands are similar to 
the COBE/DIRBE 1.25 and 2.2 ^m fihers, however the 
DIRBE filters are 0.05 fj,m (15%) wider at the long wave- 
length end. Similarly, the COBE/DIRBE 3.5 fim filter 
is slightly bluer and about 25% wider than the Spitzer 
3.6 fim filter. These slight filter miss-matches may in- 
crease the uncertainties on the measured AGB and RSG 
flux fractions by modest amounts. There are no filter 
mis-matches for the 3.6 — 24 /im fiux ratios, as all mea- 
surements are made from the same Spitzer images. 

2.4. Uncertainties 

In the following section we report the fractional flux 
contribution of TP-AGB and RHeB stars to the Magel- 
lanic Clouds as given by: 



Flux Fraction : 



Stellar Flux — Foreground Stellar Flux 
Integrated Galaxy Flux 

The uncertainties for these measurements are con- 
structed from the following components. First, we as- 
sume the published uncertainties for the integrated fluxes 
of the Magellanic Clouds at each wavelength. Second, 
we compile the photometric uncertainty in the integrated 
stellar flux from the photometric uncertainty of each star 
added in quadrature. Third, we estimate the uncertainty 
on the foreground correction by taking the Poisson uncer- 
tainty on the numbers of foreground stars in each clas- 
siflcation multiplied by the typical flux of those stars. 
We propagate these uncertainties through Equation 1. 
Fourth, we include an estimate of the uncertainty intro- 
duced by AGB stars that are fainter than the TRGB as 
described in section 2.1. 

Finally, we have quantified how much of the TP-AGB 
and RSG population we have missed due to limited spa- 
tial coverage within the apertures we used to measure 
each galaxy's integrated flux. In both galaxies, the TP- 
AGB and RSG stars follow a spherical distribution cen- 
tered on the Bars (e.g., Blum et al. 2006; Boyer et al. 
2011). For the SMC, the SAGE footprint includes an ad- 



ditional strip of coverage that extends >6 degrees along 
the SMC Tail, and for the LMC, the SAGE coverage 
stretches along the Disk, >4.5 degrees from the Bar's 
center. We can use this extended coverage to measure 
the radial profile of AGB stars in the Bars by counting 
the number of AGB stars in concentric ellipses centered 
on the Bars. We find that the AGB population is rather 
tightly concentrated, dropping to very small numbers by 
about 2 degrees (along the minor axes) from the popu- 
lation centers. By extrapolating these radial profiles in 
all directions, we estimate that we have missed about 
500 SMC AGB stars and 900 LMC AGB stars (10 SMC 
RSGs and 20 LMC RSGs) by restricting our apertures 
to 2.25 and 3.75 degrees for the SMC and LMC, respec- 
tively. Assuming that the missed stars have a similar flux 
distribution as the measured TP-AGB/RSG population, 
this amounts to an extra <1% uncertainty on the inte- 
grated AGB (and RSG) flux, which we have included in 
the upper uncertainties in Table 1 and Figure 2. 

3. RESULTS: THE MULTI- WAVELENGTH CONTRIBUTION 
OF TP-AGB AND RSG STARS TO THE INFRARED 
FLUX OF THE MAGELLANIC CLOUDS 

Figure [2] shows the fraction of the integrated galaxy 
light produced by TP-AGB (open symbols) and RSG 
(plus symbols) stars in the SMC and LMC as a function 
of wavelength from 1-24 fj,m. These results are also pre- 
sented in Table 1. At 3.6 /im, these rare stars contribute 
over 30% of the SMC flux even though they represent a 
negligible fraction of the stellar mass of the galaxy (num- 
bers of stars of each type are also shown in Table 1). In 
the LMC, they contribute more than a quarter of the 
3.6 /im light. At 24 /im, the stellar contribution to the 
total luminosity is <3% for both galaxies. 

The flux contribution from RSG stars, which are bluer 
than AGB stars, peaks near 2 /im. In the SMC, they 
are responsible for over 10% of the K-hand light. In the 
LMC, the fraction is lower, 7%. Part of this difference is 
no doubt due to the detailed star formation histories of 
these two galaxies over the last 0.1 Gyr, with the SMC 
having a higher specific star formation rate over that 
period fe.g. [Harris fc Zaritskv.2004 ,200a) . 

The flux contribution from TP-AGB stars peaks in the 
Spitzer IRAC bands between 3 and 4 /xm, reaching 15 — 
20% for both the LMC and SMC. The similarity of the 
results for the two galaxies suggests that over the longer 
time scale of 0.1 — 2 Gyr, these galaxies have had similar 
star formation histories, as has been in f erred by studies 
of their CMDs ([Harris fc Zaritskvl[200l [20091) . ' 

Figure [2] also shows the TP-AGB populations sub- 
divided into the 3 categories previously discussed: O- 
AGB, C-AGB, and x-AGB. Despite the 'fact that the O- 
AGB population outnumbers the C-AGB population by 
a factor of 1.8 in the SMC and 2.4 in the LMC, the con- 
tribution to the global fiux of the 0-AGB and C-AGB 
populations are similar at A > 3 /xm, with each con- 
tributing roughly equal parts to the total TP-AGB fiux. 
However, the contribution from 0-AGB stars peaks at 
shorter wavelengths (near 2.2 /im) than the C-AGB star 
contribution (near 3.6 /im). In contrast, the x-AGB stars 
account for less than 1% of the global flux at A < 2 /xm, 
but exceed the flux contribution of the other two TP- 
AGB types at A > 4 /xm. This is true even though the 
x-AGB population is tiny by numbers, comprising <4% 
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Fig. 2. — The fractional contribution from TP-AGB stars (open symbols) and RSG stars (plus symbols) to the integrated near- and 
mid-IR fluxes of the SMC and LMC (top and middle panels). The contributions divided by sub-species are given in the middle panels, 
while the total contributions from all TP-AGB stars (circles) and TP-AGB + RSG (asterix) arc shown in the t op panels. The bottom 
panel shows the integrated SEDs of the LMC and SMC and a model of the dust emission from Bot ct aT] l l20Tot) . While the stellar flux 
declines with wavelength (dotted line), the dust emission increases. The flux contribution of TP-AGB stars peaks at ~ 3 — 4 fim. The 
contribution from RSG stars peaks at shorter wavelengths (roughly 2.2 /im). Combined, these evolved stars contribute 32% of the SMC 
light at 3 — 4 fim and 19% at 8 fim. The flux contributions of the TP-AGB and RSG at 24 fim are small, < 3%, in both galaxies. 



of the total number of TP-AGB stars in both galaxies. 

4. DISCUSSION 

We present, for the first time, the near- through mid- 
IR flux contributions of complete samples of TP-AGB 
and RSG stars to the integrated luminosities of the 
Magellanic Clouds. We find that these stars can con- 
tribute non-negligible flux fractions from 1 to 8 fj,m (e.g., 
15 — 30%), and smaller fractions at 24 /im(e.g., < 3%). 
Below we compare these results to previous work and 
discuss the implications for estimating stellar masses and 
star formation rates of galaxies. 

4.1. The Importance of TP-AGB and RSG Stars at 
Near-IR Wavelengths 

Rest-frame near-IR fluxes typically anchor th e stel- 
lar mass estimates of galaxies (e.g., [M arasto n et al 
2006l:lConrov et al.l2009tlZibetti et al.i2009tiEskew et al 



20121 ). At these wavelengths, dust obscuration is sig- 
nificantly reduced compared with optical light, and a 
galaxy's near-IR luminosity is much less affected by the 
short-lived O- and B-type stars produced by current star 
formation. Instead, the near-IR is typically thought to 
be dominated by older, easily modeled RGB stars. 

However, we have shown that TP-AGB and massive 
RSG stars can also account for large fractions of the 



near-IR light (e.g. > 30% at 3.6 ^im for the SMC), 
even though they contribute essentially nothing to the 
total stellar masses of galaxies. A galaxy with the same 
3.6 pLTn luminosity as the SMC but comprised of entirely 
old stellar populations (e.g., > 3 Gyr) must be nearly 
twice as massive as the SMC to make up for the missing 
stellar luminosity provided by the TP-AGB and RSG at 
that wavelength. Thus, a galaxy's flux at 3-4 /im (and 
hence its M/L ratio) is highly sensitive to the SFH. 

The large TP-AGB and RSG near-IR flux contribu- 
tions we measure for the Ma gellanic Clouds fit well w ith 
the contributions reported in lMelbourne et all ()2012D for 
a sample of 23 nearby dwarfs and small spirals, which 
ranged from 5 — 40% at 1.6 /zm. Simil arly, our result s 
are well matched to the predictions of iBruzuaJ) (|2011| ). 
which discussed the contribution of AGB stars to the 
i^-band flux of their stellar population synthesis mod- 
els. Their results for a galaxy with Z — 0.008 and a 
declining SFR with an e-folding time of 1 Gyr are con- 
sistent with what we see in the Magellanic Clouds. These 
large flux contributions suggest that galaxies with even 
modest populations of TP-AGB and RSG stars will have 
very different near-IR M/L ratios than galaxies lacking 
these populations, making near-IR based stellar mass es- 
timates difficult. 

However, some authors find significantly smaller near- 
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Fig. 3. — Histograms of the [3.6] - [4.5] (top) and K3 - [8] (bottom) colors of the AGB (black) and RSG (yellow) stars in the LMC. Also 
shown are va rious color selections used in the literature to select AGB populations. The top panel shows the [3.6] — [4.5] color range used by 
IMeidt et al.l Ij2012b ) (shaded region) to isolate AGB and RSG stars, which roughly matches the mean selection used in lGerke fc KochaneM 
II2OI2I) . Both of these selections are significantly redder than the mean AGB and RSG colors in the LMC (dotted line; and the SMC, see 
Boyer et al. 2011). The lower panel shows the typical carbon star color used in the Kelson fc Holden (201 0) models (dashed line), which 
is significantly redder than the mean color for the x-AGB stars (green dash-dot line) and over 4 magnitudes redder than the mean carbon 
star color (red dash-dot-dot-dot line). 



IR flux contributions from TP- AGB stars. For instance, 
iGerke fc Koclianck (2012') calculate a smaller 3.6 /im flux 
contribution from the TP- AGB in the IVIagellanic clouds, 
only 5% of the integrated light. They selected AGB 
stars using the [3.6] — [4.5] color, however, they selected a 
CMD region that is comprised predominantly of x-AGB 
stars. Thus they miss the bulk of the TP- AGB alto- 
gether. The top panel of Figure [3] shows the distribution 
of [3.6] - [4.5] color s of T P- AGB stars in the LMC. The 
iGerke fc KocharieS (|2012l ) sample has a mean color sim- 
ilar to the x-AGB alone (green dot-dashed line), which 
are signiflcantly redder than the overall sample (dotted 
line). If we compare only the x-AGB stars in our sam- 
ple to their result, we flnd a similarly small contribution 
of ^ 3%. Including the full TP-AGB sample results in 
a much higher flux contribution at 3.6 /im, as discussed 



a bove. 

'IMeidt et al.' f2012a") shows a similar result to 
Gcrkc fc Kochan ck (2012_i), only they used unresolved 
Spitzer 3.6 and 4.5 /xm observations of galaxies in the 
Spitzer Survey o f Stellar Structure in Galaxies (S^G 
iSheth et al.ll201(il ). Their work spatially separates (on a 
pixel-by-pixel basis) older stellar populations from warm 
dust, PAH emission, and evolved stars, based on the col- 
ors of regions in the Spitzer data. Thus they are able 
to nicely constrain the warm dust and PAH emission to 
roughly 5 — 15% of the 3.6 /zm flux in a sample of 6 
nearby galaxies. Likewise they find that pixels domi- 
nated by evolved AGB and RSG stars account for only 
another roughly 5 % of the 3.6 Atm light. However, as 
with the Gerke fc Kochanekl (|2012f) result, their assumed 
AGB color of [3.6] - [4.5] = 0.24 - 0.76 (grey shaded re- 
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gion of Figure ^ is significantly redder than the bulk of 
the TP-AGB stars in the LMC which have a mean color 
of [3.6] - [4.5] = 0.0 (red dot-dot-dot-dashed line). In 
fact, the mean [3.6] - [4.5] color for TP-AGB and RSG 
in the L MC is [3.6] - [ 4.5] = -0.07 (dotted line), exactly 
the color iMeidt et all (|2012aD assign to old sta rs. Unfor- 
tunately, that means that iMeidt et al.l ()2012af) has only 
selected the reddest TP-AGB stars and, as a result, they 
are likely under-estimating the contribution of TP-AGB 
stars at 3.6 /im. 

The results we have presented here suggest that the 
near-IR may not be optimized for stellar mass estimates 
of galaxies. The 3.6 fim band, in particular, is the most 
problematic for the effects of short-lived TP-AGB and 
RSG stars, which can impose a stochastic variability 
on the near-IR M/L ratio of stellar populations (e.g., 
[Melbourne et al.ll2012[ ). However, it should be possible 
to minimize the effects of the TP-AGB and RSG stars on 
estimates of the near-IR M/L ratios of galaxies by using 
bluer wavebands. The TP-AGB and RSGs flux contri- 
bution at 1 /im is roughly half that of the 3.6-/im flux 
contribution. Of course, bluer wavelengths will be more 
impacted by extinction from dust and emission from hot 
O and B stars, but the /- or J-band may prove to be 
the ideal compromise. To test this hypothesis we com- 
pare the point-source flux above the TRGB to the total 
point-source flux (which is measurable to 2.5 mags below 
the TRGB) as a function of wavelength from 0.5—3.6 /im. 
Figure m displays the results, and shows that indeed there 
is a minimum flux contribution from luminous short-lived 
stars as a function of wavelength. The minimum occurs 
between 0.8 and 1 /im. At shorter wavelengths, massive 
main sequence O and B stars contribute large fractions 
of light. At longer wavelengths, the TP-AGB and RSG 
stars contribute large fractions. This plot suggests that 
rest- frame /— through J— band will indeed provide M/L 
ratios that are the most stable against stellar population 
history and the influence of these short-lived luminous 
stars. 

4.2. The Importance of TP-AGB and RSG Stars at 
Mid-IR Wavelengths 

Because mid-IR light is produced by warm dust and 
PAH emission in star forming regions, mid-IR luminosi- 
ties of galaxies are often used as proxies for their cur- 
rent star formation rates. However, we have shown that 
the TP-AGB and RSG stars can also contribute to the 
mid-IR fluxes of galaxies, with the bulk of the mid-IR 
contribution coming from the rare but extremely dusty 
x-AGB stars. In the SMC, the contribution of TP-AGIB 
and RSG stars at 8 /im is wl9%! If these x-AGB stars 
are neglected in the accounting of the mid-IR light, an 
SMC SFR based on 8 /im light alone will be biased high. 

This result has impacts on the interpretation of thou- 
sands of mid-IR based star formation rates, for example 
at z = 2, Spitzer 24 /im imaging traces rest-frame 8 /im, 
where the TP-AGB is potentially an important flux con- 
tributor. However, the size of the effect may be small 
in practice. For instance, in the LMC, the TP-AGB + 
RSG 8 /im contribution is < 5%. Compared with the 
SMC, the LMC has a significantly lar ger contribution 
from warm dust and 8 /im PAH emission ( Bot et aO2010l 
and Figure[2]here). Samples of the extremely dusty z = 2 
galaxies detected by Spitzer have been shown to have 




1 

Wavelength {/^m) 



10 



Fig. 4. — The contribution of luminous stars (with fluxes > than 
the TRGB) to the total point-source flux (measured to 2.5 mags 
below the the TRGB) as a function of wavelength. At the blue end 
the O and B stars contribute significantly to the point-source flux. 
At the red end, RSG and TP-AGB stars contribute significantly to 
the total. As all of these stars are short lived, they will result in a 
rapidly changing M/L ratio at these wavelengths. The impact of 
luminous short-lived stars on the stellar mass estimates of galaxies 
can be minimized by using the M/L ratio at 0.8 — 1.0 fim, where 
this plot shows a minimum. 



large PAH flux contributions fe.g. iDesai et al.ll2009( ). so 
the correction for TP-AGB stars may be small in the 
high-z luminous infrared galaxies. However, it is possi- 
ble that TP-AGB stars artificially elevate the the rest- 
frame 8 /im estimated SFRs of galaxies, especially for 
lower-metallicity systems like the SMC where (I) PAH 
emission is limited, (2) C-AGB stars are more common, 
and (3) star formation rates are likely smaller. 

We are not the first to point out that the TP-AGB 
may be importan t at mid-IR wavele ngths. For in s tance , 
iKelson fc HoldCTl UoiO), and Chis ari fc KelsOTll (poll 
recently proposed TP-AGB stars as the primary source 
of mid-IR flux in galaxies, and may account for the rise 
of luminous i nfrare d galaxies with redshift. Similarly, 
iVerlev et al.l ()2009D suggested dusty shells surrounding 
AGB stars could account for the bulk of diffuse 24 /im 
emission in the outskirts of M33. While our results show 
the importance of TP-AGB stars at mid-IR wavelengths, 
these earlier works likely overemphasize the TP-AGB 
contribution, especially at 24 /tm, where we show that 
TP-AGB stars contribute less than 3% of the Magellanic 
Cloud flux. 

The large mid-IR TP-AGB contributions predicted 
by these previous works (e.g. > 50% of the 8 and 
24 /im flux) likely arose from their adopted mid-IR 
colors and fluxes of C-AGB stars . IKelson fc Holdeiil 
((2010h and lChisari fc KelsonI (l20T2h assume that C-rich 
AGB stars have an average K — [8.8] color of 5.66 ± 
1.16 mag, derive d from Galactic c a rbon st ars observed 
in the mid-IR bv iGuandalini et al.l (|2006l ). IVerlev et afl 
l|200l assume an average 24 /tm flux derived from 
carbon st ars observed with Spitze r in the Magellanic 
Clouds bvlGroenewegen et al. (I2007D. However, b o th the 
IGuandalini et al.l ()2006f ) and ICroenewegen et al.l (|2007D 
samples are heavily biased towards the reddest x-AGB 



8 



stars and completely exclude the bulk of the C-AGB pop- 
ulation, which is much bluer. This issue is demonstrated 
in the bottom panel of Figure |3] which plots histograms 
of the Ks - [8] color of the TP-AGB stars in the LMC. 
The color adopted by lKelson fc H olden' (2010) models is 
shown as a dashed line and is significantly redder than 
the bulk of the LMC C-AGB stars. When including the 
complete samples of C-rich AGB stars (C-AGB plus x- 
AGB, red and green distributions respectively) from the 
SAGE survey, we find {Kg — [8])x+c = 1-38 mag in the 
SMC and 1.31 mag in the LMC (dot-dot-dot-dashed line 
in Figure [3]). This is over 4 magnitu des bluer than the 
Kelso n & Holden models (also used in iChisari fc KelsonI 
I2012D . A similarly sized discrepancy is obtained when 
comparing the SAGE AGB sample to the adopted 24 fim 
flux for AGB stars in Verley et al. These revised colors 
would result in a much smaller contribution of TP-AGB 
stars at mid-IR wavelengths, as is seen Figure [2] for the 
Magellanic Clouds. 

While some x-AGB stars in the SAGE Magellanic 
Cloud sample do show colors as red as those used in the 
Kelson & Holden models, together they have an average 
color that is still 2 magnitudes bluer than the models, 
{Ks - [8])x, SMC = 3.25 mag and {Ks - [8])^, lmc = 
3.41 mag (dot-dashed line in top panel of Figure In 
addition, these x-AGB stars contribute almost nothing to 
the K-band light of the Magellanic clouds (1% or less). 
The bluer 0-AGB and C-AGB stars produce roughly 20 
times the ii'-band light of the x-AGB stars. The Kelson 
& Holden models rely on the X-band flux to estimate the 
AGB contribution at longer wavelengths, and they as- 
sume that the x-ABG stars are contributing nearly 50% 
of the fC-band light. Because of these assumptions, they 
are likely significantly over-estimating the AGB flux con- 
tribution at longer wavelengths. 



5. CONCLUSIONS 

The contribution of TP-AGB and RSG stars to the 
integrated fluxes of the Magellanic Clouds, varies with 
wavelength. We find that their combined contribution 
peaks between 3 — 4 fj,m, where they contribute 32% of 
the SMC flux and - 25% of the LMC flux. At 1.2 ^m, 
the contributions are smaller, ^ 18% and ^ 15% for the 



SMC and LMC. At longer wavelengths the contributions 
also decline, but still constitute 19% of the 8 /im light of 
the SMC (only 4% for the LMC where warm dust and 
PAH emission are stronger). The bulk of this 8 /im stellar 
flux is from the most dust-enshrouded AGB stars, the 
so-called x-AGB stars, which are predominantly carbon 
stars near the end of the AGB phase. At 24 fim the AGB 
star contribution to the total flux is less than 3% for each 
galaxy. 

The flux contributions of the TP-AGB and RSG are 
significant at 3.6 /im. Thus, these stars have an im- 
pact on the M/L ratios of the Magellanic Clouds at 
near-IR wavelengths. Because these stars are short-lived 
we can expect that the near-IR M/L ratios will also 
vary significantly over time. To minimize the impact of 
short lived TP-AGB and RSG stars on stellar mass es- 
timates of galaxies, shorter wavelength observations are 
preferred. However, at optical wavelengths, short-lived, 
super-luminous O and B stars play a similar role as the 
AGB and RSG stars in the near-IR. We find that the op- 
timal wavelengths for minimizing the impact of luminous 
short-lived stars on the M/L ratio of galaxies is between 
0.8 and 1 /im. 

The non-negligible fiux from the TP-AGB at 8 /im sug- 
gests that star formation rates based on rest-frame 8 /tm 
flux (e.g. Spitzer 24 /tm imaging at z = 2) may be bi- 
assed slightly high. However, this effect is most impor- 
tant for systems with little PAH emission (like the SMC). 
For systems with significant warm dust and PAH emis- 
sion, the TP-AGB contribution at 8 jiva appears to be 
small, thus for the LMC the contribution is <4%. 

Special thanks to Karl Gordon for providing integrated 
fiuxes of the LMC in the Spitzer bands, to Caroline Bot 
for sharing the SED and model fits of the Magellanic 
Clouds. Thanks to Gustavo Bruzual for sharing updated 
SSP models of the Magellanic Clouds. We also want to 
thank Barry Madore, Mark Seibcrt, and Wendy Freed- 
man for inspiring efforts to study the impact of AGB 
stars on mid-infrared galaxy fiuxes. This work is based 
in part on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology under a contract 
with NASA. 
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TABLE 1 

Integrated Flux Density of the SMC & LMC from 1 — 24 /im 
AND Percentage of that Flux Contributed by TP-AGB and RSG Stars 



Type 


Stars^ 


Fgnd"^ 


1.2^tm<^ 


2.2^tm<^ 


3.5Mm<i 


3.6/jm'^ 4.5/xm° 


5.8Atm° 


8Atm<= 


24^m^ 




SMC 












Integrated Flux [Jy] 














670 ± 80* 


525 ± 65'' 


280 ± 40'' 


300 ± 20s 220 ± lOS 
Fractional Contribution [%] 


220 ± 108 


200 ± lOS 


350 ± lOS 


O-AGB 
C-AGB 
x-AGB 
RSG 


3353 
1597 
309 
2340 


4 

< 1 

< 1 
598 


5 7+1-1 
-0.7 

3 O+O-" 


04-1-1.5 
''•^-i.o 
6 6+° '' 

10.911:1 


8.0+li 
9 1+^-* 

4.8«-; 

w.otii 


7 t;+i.l 6 0+°-^ 
'-3-0.5 "-^-0.3 
+0.8 71+0.6 

4 ^+°-^ 6 ^+°-^ 
*-3_0.3 "-"^-0.3 
04+0.6 7 6+0-3 
^-^-0.6 '-"-0.3 


4 ^+0.6 
^-"-0.2 
c 0+0.4 

0.0_Q 2 

"-="-0.3 

r r+0.3 
■^-"-0.3 


1+0.4 
■'-1-0.2 

^■''-0.2 

'-"-0.4 

.1 7+0.2 
''■ ' -0.2 


14+0-01 
"-1^-0.01 

31+0-02 
""'1-0.01 

1 30+0-04 
l-''"-0.04 
Q 00+0.01 
""'''-0.01 


Total 


7599 


602 


i7.8ti:«3 


27 1 +'^- ' 


31.9ltt 


29.81i;q 27.2loj 


90 + 2. u 
^^-■^-0.5 


19-3^0:5 


9 -I+U.2 
^•1-0.1 


LMC 












Integrated Flux [Jy] 














4520 ± 650f 


3770 ± 540' 


2190 ± 300' 


2080 ± 100^ 1350 ± 70^ 
Fractional Contribution [%] 


1850 ± 200^ 


5980 ± 300^ 


7640 ± 200^ 


O-AGB 
C-AGB 
x-AGB 
RSG 


15635 
5628 
989 
4361 


< 1 

< 1 

< 1 
13 


6.8l?:« 

^-0-0.4 

n 1+°-^ 

'-'■1-0.1 

rr ^ + 0.8 


9.9+l\ 

6 n+°-^ 

7+°-^ 
'J-'-o.i 

D.y_Q 7 


« s+3-1 
°-°-0.8 

7 6+°-^ 
'•"-0.7 
2 q+0.3 
^•^-0.3 

6 n+°-^ 
"•'-'-0.5 


Q + 3.2 Q + 3.1 

^-■^-0.5 "-■'-0.4 
1+0.5 7 0+0.4 
"-1-0.4 '-'-'-0.4 
r,+0.2 c 9+0.3 
"-'^-O.l ■-'-^-0.3 
c 0+0.3 c c:+0.3 
"-'-'-0.3 ■-'-"-'-0.3 


7+0.4 
' -0.4 
4 0+0.5 

Q+0.3 
'^-'-'-0.3 


9+0-3 

1 0+0-1 
i-'^-o.i 


n+o-22 
"-1I-0.01 

n n+0-01 
"-11-0.01 

3'S+O-Ol 
"-■'O-O.Ol 

n 16+0-01 

"-1"-0.01 


Total 


26613 


13 




23.5li:4 


25. 2l^;^ 


26.6lo7 26.3lo7 


15 6+^" 
10-0-0.9 


3-9lii:t 


"•'-0.1 



^ Number of stars of given typo detected at 3.6 /im; the number varies slightly -with -wavelength 
^ Number of Milky Way foreground stars as estimated by TRILEGAL tGirardi ct al..,2005h 

Integrated Flux from COBE/DIRBE (at 1.25 and 2.2 fim), Stellar Flux from 2MASS (at 1.2 and 2.2 /jm) 
'I Integrated Flux from COBE/DIRBE (at 3.5 /im), Stellar Flux from Spitzer (at 3.6 iJ,m) 
^ Integrated Flux from Spitzer, Stellar Flux from Spitzer 
1 From Israel ot al. (20 10) 
S From lGordon et al.l 1201 ll ) 
1' From Gordon et al. (in preparation) 



